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Abstract 

We study dynamics of two coupled periodically driven oscillators. Im- 
portant example of such a system is a dynamic vibration absorber which 
consists of a small mass attached to the primary vibrating system of a 
large mass. 

Periodic solutions of the approximate effective equation are determined 
within the Krylov-Bogoliubov-Mitropolsky approach to get the amplitude 
profiles A{^1). 

Dependence of the amplitude A of nonlinear resonances on the fre- 
quency n is much more complicated than in the case of one Duffing os- 
cillator and hence new nonlinear phenomena are possible. In the present 
paper we study metamorphoses of the function A (57) induced by changes 
of the control parameters. 

1 Introduction 

Coupled oscillators play important role in many scientific fields, e.g. biology, 
electronics, and mechanics, see [1] [H |3] and references therein. In this paper we 
analyse two coupled oscillators, one of which is driven by an external periodic 
force. Important example of such system is a dynamic vibration absorber which 
consists of a mass m2, attached to the primary vibrating system of mass mi 
[H [S] . Equations describing dynamics of such system are of form: 



mixi - Vi (xi) - Ri (xi) + V2 (±2 - ii) + i?2 {^2 - xi) = / cos (ujt) 
1112X2 - V2 {X2 - Xi) ~ R2 {X2 - Xi) ^ 



(1) 
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where Vi, Ri and V2, R2 represent (nonlinear) force of internal friction and 
(nonlinear) elastic restoring force for mass mi and mass 7712, respectively. In 
the present paper we do not assume that the ratio 1712 /rni is small. 
In the present paper we shall consider a simplified model: 

Ri {xi) = -aixi, Vi (ii) = -vixi. (2) 

Dynamics of coupled periodically driven oscillators is very complicated [TJ 
m [3]. We simplified the set equations ([T]), ([2]) by reducing it to the prob- 
lem of motion of two independent oscillators. More exactly, we derived the 
exact fourth-order nonlinear equation for internal motion as well as approxi- 
mate second-order effective equation in 0. Moreover, applying the Krylov- 
Bogoliubov-Mitropolsky method to these equations we have computed the cor- 
responding nonlinear resonances (cf. [6] for the case of the effective equation). 
Dependence of the amplitude A of nonlinear resonances on the frequency uj is 
much more complicated than in the case of Duffing oscillator and hence new non- 
linear phenomena are possible. In the present paper we study metamorphoses 
of the function A (cj) induced by changes of the control parameters. 

The paper is organized as follows. In the next Section derivation of the exact 
4th-order equation for the internal motion and approximate 2nd-order effective 
equations in non-dimensional form are presented. In Section 3 metamorphoses 
of amplitude profiles determined within the Krylov-Bogoliubov-Mitropolsky ap- 
proach for the approximate 2nd-order effective equation are studied and the case 
of the standard Duffing equation is presented as well. More exactly, theory of 
algebraic curves is used to compute singular points on effective equation am- 
plitude profiles - metamorphoses of amplitude profiles occur in neighbourhoods 
of such points. In Section 4 examples of analytical and numerical computa- 
tions are presented for the effective equation. Our results are summarized and 
perspectives of further studies are described in the last Section. 



2 Exact equation for internal motion and its ap- 
proximations 

In new variables, x = xi, y = X2 ^ xi, equations ([T]), ^ can be written as: 

mx + i/x + ax + Ve (y) + Re {y) = f cos (ujt) 
me {i + y)- Ve (y) - Re (y) = 

where m = mi, me = m2, v = i^i, a = ai, Ve = V2, Re ^ R2- 

Adding equations ^ we obtain important relation between variables x and 

y- 

Mx + vix + aix + mejj = f cos (ujt) , (4) 

where M = m + me . 
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We can eliminate variable x in ^ to obtain the following exact equation for 
relative motion: 

(^^ ^ + + ") (mJ/ - K (y) - Re iv)) + Ame {vf^ + a) y = Fcos (c^i) , (5) 

where F = rUeUp' f , ^J^ = mrric/M and A = nie/M is a nondimensional parameter. 
Equations ([S]), (|1]) are equivalent to the initial equations (P), © [B]. 
In the present work we assume: 

Re (y) = tteV - '^eV^ , Ve (jf) = -I^ejf. (6) 

We thus get: 

f^^^ + + f + >^e^-aey + Jey') +Am, (,yA + a) ^ ^ F cos (..i) , 



. 

We shall write Eq. ([7]) in nondimensional form. Introducing nondimensional 
time T and rescaling variable y: 

T = tu!, z = yw — , (8) 

V 

where: 



(9) 



we get the exact equation for motion of mass mg: 

{£.+H^^+a)[^ + hj^^-z + z') +n {Hj, + a) ^ G^17^ cos {^r) , 

(10) 

where nondimensional constants are given by: 

h=^, H=^, n = ^, G = -^^f^f, K^^, a = (11) 

We shall consider hierarchy of approximate equations arising from (jlOp . For 
small K, if, a we can reject the second term on the left in (|10p to obtain the 
approximate equation: 

(i^ + + a) + /if - z + ^ jQ' cos m , (7 ^ G^) 

(12) 

which can be integrated partly to yield the effective equation: 

-r^ + h- z + z"^ = = cos(f)r + (5), (13) 

where transient states has been omitted jSj. And finally, for = 0, a = we 
get the Duffing equation: 

^ + h^-z + z^ = -^Gco8{nT + S). (14) 
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3 Metamorphoses of the ampHtude profiles 

We applied the Krylov-Bogoliubov-Mitropolsky (KBM) perturbation approach 
[3 [8] to the effective equation (fT3)) obtaining for the 1 : 1 resonance the foUowing 
amphtude profile ^6j: 

^eff - =■ (15) 

Now, for iJ = 0, a = 0, we obtain the amplitude profile for the Duffing equation 
" - ^ (16) 



h^n^ + (1 + r!2 - ^Al) 



It is well known that dependence of the function Ad, cf. p^ . on control 
parameters 7, h is rather simple. On the other hand, dependence of the ampli- 
tude profile Aeff (ft) on control parameters 7, h, a, H is more complicated and 
thus Aeff (17) can describe new nonlinear phenomena. In the next Section we 
shall study possible metamorphoses oi Ad, A^ff induced by changes of control 
parameters, the more complicated case of the 4th-order exact equation (ITU|) will 
be treated elsewhere. 

Equations (jl6p . (I15|) define the corresponding amplitude profiles implicitly. 
Such amplitude profiles can be classified as planar algebraic curves. Firstly, we 
shall collect useful theorems on implicit functions which will be used below. 

Let us write equations ([T5|) . (fTH]) as Le {Ye,X) = and Ld {Yd,X) = 0, 
respectively, where X = il^, Y = A^: 

(h^X+{l+X~lY,y^(^iX-af + H^X^Y,-j^X^ = 0, (17) 
h^X+il + X^lYof^Yo-j^ = 0. (18) 

It follows from general theory of implicit functions [11] [12] that conditions 
for critical points of Y (X) read: 

Moreover, critical points of the inverse function X (Y) are given by: 

dL(Y,X) /dL(Y,X) \ 

L(Y,X)^Q, QY \ dX ^ ^)- 

It may happen that in some points {Xo,Yq) we have: 

L{Y,X) = Q, ^^^^=0, L{Y,X)=Q. (21) 

Such points are referred to as singular points of algebraic curve L {Y, X) = 
because they are in some sense exceptional. 



4 



3.1 The case of the DufRng equation 

Singular points of the algebraic curve defined by are given by: 



dX 






(22) 




(23) 


0, 


(24) 


0, 


(25) 


0, 


(26) 




(27) 



dY 

The set of equations (IT51) . (P^ . ([25]) can be written as: 

h^XY + Y {I + X - ^Yf ~ -i'' 
h^Y + 2Y + 2YX - 
h'^X + 1 + 2X ^3Y + X^ - 3YX + ^Y^ 

where X, Y are positive. 
General solution reads: 

+ |r (y ^ 0) 

It follows that Y > 0, X < and 7^ < and thus the system of equations 
([M]) . (|26p has no acceptable solutions since we assume that h, 7 are real 

and X, Y are non-negative. 

3.2 The case of the effective equation 

Singular points of the algebraic curve defined by (|17|) are given by equations: 

w = <^«' 

It follows from (IT7)) and ((29|) that either of equations must hold: 

16X2 + le/i^X + 32X + 16 - 48rX - 48r + 27^2 = 0, (30a) 

{X-af + H^X = 0, (30b) 

where X, Y are positive. 

Let us start with Eq. (|30b| . In this case we obtain from dTT]), ([Ml) and (|30b| 
the following rather special solution: 

X = a, i7 = 0, 7 = (y, /i, a - arbitrary) . (31) 

Let us now consider more general Ea. pOa|) . We can treat X as arbitrary. 
Then we obtain two solutions for Y: 



Y = I + fx - I Vl + 2X + X^- 3h^X, (32) 

Y = f + fx + Wl + 2X + X^- iK^X, (33) 
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where the inequaUty 

1 + 2X + X^ - 3h^X > 0, (34) 
must hold. This means that for a chosen value of X the parameter h must obey 



3X 



Solving equations dTT]), ([Ml), and ^ or ^ we get: 

Y — 8j8\^_i_4j 

a = Zi- 
H 



IX±p{X), 

X 

-3- 



(h'^X+X^+2X + l)^h 



(ft2x+X2+2X + l)3/i 



where U and Zi, Z2 are given by: 

u{x) = yr 



2X + X^ - 3h^X 
Zi = Vwi iX)±W2 {X)U{X), 
Z2 = VW3 iX)±W4 iX)UiX), 



with 



wi {X) 




+ 05^^ - 


f a4X** + a3X^ + 02^2 + aiX + oq 




0-6 


= 16/i2, 








05 


= ASh"^ - 


f 96/1^, 






04 


= 48/1^ - 


f 192/1^^ H 


-240/42 + 672, 


< 


03 


= 16h^ - 


f 96/1^ + 


288/1^ + (320 + 10272) + 2472, 
- (I6272 + 240) /i2 +3672, 
+ (I872 + 96) /i2 +2472, 
+ 672, 




02 


= 48h^ - 


f 192/1^^ H 




ai 


= (5472 


+ 48) h'^ 






= (16- 


4272) 



u;2 {X) = &3X3 + 62^2 + biX + 60, 
63 = 672, 

62 = -5I72/12 + I872, 
bi = -972/1^ _ 1272/12 _^ 13^2^ 

60 = 3972/12 + 672, 



W3 (X) = crX^ + ce^e 
-32/i2, 



- 05^5 + CiX^ + CgX^ + C2X-^ 



(35) 

(36) 
(37) 

(38) 

(39) 
(40) 
(41) 

(42) 



(43) 

(44) 
(45) 

ciX + co, (46) 



C7 
C6 
C5 
C4 

C3 

C2 
Cl 
Co 



-96h^ 

-32/i« 
-640- 



- 192/i2, 

- 384/1* _ 
- 192/1^ - 

4272) /l2 



480/i2 + 32Zi/i, 
576/1^^ + 64Zi/i3+ 
+ 128^1/1 + 672, 



: -96/1^ + 32Zih^ + (-384 + 5472) 
USZih^ + (-480 + I872) /i2 + 192Zi/i + 2472, 
= I6272/12 + 64Zi/i3 + 128Zi/i - 192/i2 + 3672 - 96/i*, 
= (10272 - 32) /i2 + 32Zi/i + 2472, 
= 67^ 

u;4 (X) = &o^^ + + b2X + 63. 



(47) 



(48) 
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4 Analytical and numerical computations 



It follows from solutions obtained in the preceding Section that we can control 
position of a singular point. More exactly, we choose a value of X and then 
h fulfilling inequality ([55]) can be chosen as well. Next we specify 7 and then 
y, a, H are computed from Eqs. (155)) . ([57|. psp . In this process the position 
of the singular point (X, 1") and values of control parameters H, h, 7, a are 
determined (provided that the solutions are real) . 

Bifurcation diagram for the effective equation (fT3)) for the following values 
of control parameters H = 0.04, h — 0.4, a — 0.8, 7 = 2.5 is shown in Fig. [T] 
(cf. Fig. 1 in [6]) where colours mark different initial conditions. Position of the 



-_z 
1 




-3 



- , , . ,0,5. , , ,1, . , , ,1,5, . , ,2 , , . , ,2.5, , . ,3, , , . ,3,5. CI 



Figure 1: Bifurcation diagram for the effective equation, h = 0.4, 7 = 2.5, 
a = 0.8, H = 0.04. 

1 : 1 resonance agrees well with the amplitude profile, computed for the same 
parameters, cf. Fig. [5] and discussion in [5]. 

We shall now compute coordinates of a singular point using equations , 
dSZI), dMl). At first we choose the value oi X as X = 9 {n = 3). Then we 
can select any value of h obeying inequality (1351) . We thus put h ~ 0.8 to get 
from Eq. ^ Y = 4.8466 {A = 2.2015). Next we choose 7 = 1.5 to compute 
from dSIl), dSHl) a = 9.0720, H = 0.2995. In Fig. |3] below we show amplitude 
profiles computed from Eq. ()15|) for critical parameter values {H, h, a, 7) = 
(0.2995, 0.8, 9.0720, 1.5), and for two more values oi H, H > H^r and H < 

Her- 
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Q 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 

omega 

Figure 2: Amplitude profile A{n), h = 0.4, 7 = 2.5, a = 0.8, H = 0.04. 




0.0 0,2 0.4 0.6 0.8 1.0 1,2 1.4 1.6 1.8 2.0 2.2 2.4 2.fi 2.8 3.0 3.2 3.4 3.6 3.8 4.0 

omega 

Figure 3: A{n) in the singular point and in its neighbourhood, h = 0.8, 7 = 1.5, 
a = 9.0720, and H ^ 0.27 (green), H = 0.2995 (red), H = 0.33 (blue). 

The critical (red) curve intersects itself in singular point {X, Y) = (4. 8466, 9) 
or {A, ri) = (2.2015, 3). Green curve corresponds to H — 0.27 while blue 
curve has been computed for H = 0.33 (other parameter values unchanged). 
The initial amplitude from Fig. [T]was also shown (black curve). The first bi- 
furcation diagram, cf. Fig. |4l was computed for H = 0.27 and corresponds to 
the green curve in Fig. |3l We note that the small branch of the 1 : 1 reso- 
nance is discontinuous in agreement with the amplitude profile shown in Fig. 
[3] (green curve) . The next bifurcation digram. Fig. [5l has been computed for 
critical value H = 0.3019 determined numerically from Eq. (IT^ (this differs 
slightly from the critical value H — 0.2995 determined from the KBM solution 
as described above). 
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2.6 ,2.7 ,2.B ,2.9 ,3 ,3.1 ,3.2 ,3.3 



Figure 4: Bifurcation diagram, = 0.8, 7 = 1.5, a = 9.0720, H = 0.27. 
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Figure 5: Critical diagram, h = 0.8, 7 = 1.5, a = 9.0720, F = 0.3019. 



And finally, the last bifurcation diagram was computed for H = 0.33 - and 
again the small branch of the 1 : 1 resonance is continuous. 
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2,6 ,2,7 ,2,8 ,2,9 ,3. ,3..1 ,3.2 ,3.3 fl 



Figure 6: Bifurcation diagram, h — 0.8, 7 = 1.5, a = 9.0720, H = 0.33. 

It follows from results presented in Section that ioi X — 9, h — 0.8, 7 — 1.5 
there is another singular point. Indeed, we can compute Y from another of 
equations ^ to get from Eqs. dST]), 1^ Y ^ 12.9311 (A = 3.5960), a = 
9.1213, H ^ 0.5158. 




Figure 7: A{n) in the singular point and in its neighbourhood, h — 0.8, 7 = 1.5, 
a = 9.1213, and H = 0.49 (green), H = 0.5158 (red), H = 0.54 (blue). 

In Fig. [7] amplitude profiles computed from Eq. (|15p for critical parameter 
values (i7, ft,, a, 7) = (0.5158, 0.8, 9.1213, 1.5) and for two other values of H, 
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H < Her and H > Her have been shown. Bifurcation diagrams for H = 0.49 
and H = 0.54 are shown below. 



-Z 
















1 














- 


0.5 
I 
















:0.5 
















-1 
















2.B 


,2 ..85 


,2..9 


,2 ..95 


,3. 


,3 ..05 


,3..1 


,3 ..15 ^ 



Figure 8: Bifurcation diagram, h = 0.8, 7 = 1.5, a = 9.1213, H = 0.49. 
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Figure 9: Bifurcation diagram, h = 0.8, 7 = 1.5, a = 9.0720, = 0.54. 
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5 Summary and discussion 



In this work we have studied metamorphoses of amphtude profiles for the effec- 
tive equation, describing approximately dynamics of two coupled periodically 
driven oscillators. Our analysis has been analytical although based on the ap- 
proximate KBM method. 

Theory of algebraic curves has been used to compute singular points on 
effective equation amplitude profiles. It follows from general theory that meta- 
morphoses of amplitude profiles occur in neighbourhoods of such points. In 
Section 3 we have computed analytically positions of singular points for the 
amplitude profiles A (il) determined within the Krylov-Bogoliubov-Mitropolsky 
approach for the approximate 2nd-order effective equation (fT3|) . In the first case 
the singular point corresponds to self-intersection of A (il), see Fig. |31 while in 
the second case it is a isolated point, cf. Fig. [7l 

It is interesting that the solution described in Section 3 permits control of 
position of singular point: we choose arbitrary value of variable X (X^ = fl) , 
then value of the parameter h obeying inequality (j34]) is selected. Finally the 
value of the control parameter 7 is chosen and Y, a, H are computed from 
Eqs. ((36)) . (|37| . (p8| : it should be stressed that we have not come across any 
difficulties to obtain real solutions. We hope to carry full analysis of conditions 
guaranteeing existence of real solutions in our future papers. As a by-product 
we have demonstrated that there are no singular points for A (f2) computed for 
the Duffing equation in agreement with well established numerical experience. 

We have also computed numerically bifurcation diagrams in the neighbour- 
hoods of singular points and indeed dynamics of the effective equation (|13p 
changes according to metamorphoses of the corresponding amplitude profiles. 

In our future work we are going to study singular points of the amplitude 
profiles computed for the exact equation pH]) . 
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